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ABSTRACT 

Recent results have shown that a substantial fraction of high-redshift Lyman alpha 
galaxies contain considerable amounts of dust. This implies that Lyman alpha galax- 
ies are not primordial, as has been thought in the past. However, this dust has not 
been directly detected in emission; rather it has been inferred based on extinction es- 
timates from rest-frame ultraviolet (UV) and optical observations. This can be tricky, 
as both dust and old stars redden galactic spectra at the wavelengths used to infer 
dust. Measuring dust emission directly from these galaxies is thus a more accurate way 
to estimate the total dust mass, giving us real physical information on the stellar pop- 
ulations and interstellar medium (ISM) enrichment. New generation instruments such 
as the Atacama Large Millimeter Array (ALMA) and Sub-Millimeter Array (SMA), 
should be able to detect dust emission from some of these galaxies in the sub-mm. 
Using measurements of the UV spectral slopes, we derive far-infrared flux predictions 
for of a sample of 23 z ^ 4 Lyman alpha galaxies. We find that in only a few hours, 
we can detect dust emission from 39 ± 22% of our Lyman alpha galaxies. Comparing 
these results to those found from a sample of 21 Lyman break galaxies (LBGs), we find 
that LBGs are on average 60% more likely to be detected than Lya galaxies, implying 
that they are more dusty, and thus indicating an evolutionary difference between these 
objects. These observations will provide better constraints on dust in these galaxies 
than those derived from their UV and optical fluxes alone. Undeniable proof of dust 
in these galaxies could explain the larger than expected Lya equivalent widths seen 
in many Lya galaxies today. 
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1 INTRODUCTION 

Lyman alpha (Lya) galaxies have been studied in detail 
for nearly a decade now, and we are only just beginning 
to learn about the properties of individual objects. Most 
initial narrow-band studies (e.g., Hu et al. 1998; Rhoads et 
al. 2000) observed blank fields with two filters: a narrow- 
band filter (to capture the Lya emission line), and a broad- 
band filter encompassing this narrow-band (used to measure 
a narrow-band excess). Learning the detailed properties of 
individual Lya galaxies would be difficult with these data. 
However, these early studies were able to measure the equiv- 
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alent widths (EWs) of their objects, calculated from the 
ratio of narrow-band flux to broad-band flux (Malhotra & 
Rhoads 2002). 

Some of these initial results were surprising, as in many 
cases the Lya EW distribution was higher than one might 
have expected from stellar population models, as many ob- 
jects exhibited rest-frame Lya EWs over 200 A (e.g., Ku- 
dritzki et al. (2000); Malhotra & Rhoads (2002); Dawson 
et al. (2004; 2007); Shimasaku et al. (2006); Finkelstein et 
al. (2007)). A normal stellar population (defined as having 
a Salpeter (1955) initial mass function (IMF), solar metal- 
licity, and a constant star formation rate) has a maximum 
Lya EW of 260 A at 1 Myr, falling to a constant value of 95 
A at 100 Myr. Thus, while objects with EWs > 200 A are 
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possible with a normal stellar population, they would all 
have to be extremely young. Looking at a random sample of 
Lya galaxies, it would be extraordinarily rare if they all had 
ages of only a few Myr, thus something is likely happening 
to skew the EWs higher than they should be. 

There are many possible reasons why this could be oc- 
curring, and the easiest area to examine is the assumption 
of a normal stellar population. For example, if a stellar pop- 
ulation had sub-solar metallicity, it would form stars with 
hotter stellar photospheres. These stars would then pour 
out more ionizing photons, which, after interacting with the 
local interstellar medium (ISM), would in turn create more 
Lya photons. Another possibility exists, whereby if a galaxy 
had a top-heavy initial mass function, it would create more 
massive stars than a Salpeter IMF, again resulting in more 
ionizing photons, and thus more Lya photons. Both of these 
scenarios are possible in primitive, dust-free galaxies. Strong 
Lya lines also come from the active galactic nuclei (AGN) 
components of some galaxies, but previous work has shown 
that most Lya galaxies detected via narrow-band selection 
are not AGN (Malhotra et al. 2003; Wang et al. 2004). 

Lya photons are resonantly scattered by any neutral 
hydrogen present, thus in a homogeneous ISM, they will 
have to travel much further to escape the galaxy, dramati- 
cally increasing the chance of absorption by dust (Meier & 
Terlevich 1981). However, if the dust is in an inhomogeneous 
distribution, it is possible that a galaxy could appear to have 
a strong Lya EW at later ages. 

In this scenario, Lya can escape unaffected by dust if it 
is somehow screened from interacting with the dust grains. 
This can happen if the dust resides in cold, neutral chimps 
embedded in an ionized inter-clump medium. This scenario 
was first proposed by Neufeld (1991) and examined in detail 
more recently by Hansen & Oh (2006) . In this ISM geometry, 
the Lya photons get scattered by neutral hydrogen at the 
surface of these clumps, and thus never see the dust which is 
embedded deeper inside. However, the continuum photons 
suffer no such restriction, so they pass un-obstructcd into the 
interior of the clump, with many continuum photons being 
absorbed (extincted) and scattered (resulting in reddening of 
the continuum colour). The Lya photons proceed to bounce 
around from clump to clump, eventually escaping the gala^xy. 
The net result of this is an increase in the Lya - continuum 
flux ratio over that intrinsic to the stellar population. We 
would thus detect this galaxy with a higher EW, irrespective 
of its stellar population characteristics. 

We believe that this scenario may be responsible for 
some of the high-EW objects seen in narrow-band studies, 
and in our previous work (Finkelstein et al. 2009; hereafter 
F09), we found ten objects out of a sample of 14 LAEs where 
this scenario appears to bo in effect (one of which was also 
studied in Finkelstein et al. 2008a). Two of these objects 
are evolved galaxies with large Lya EWs due to this dust 
enhancement scenario (the other eight are still young, but 
have their Lya EWs enhanced to be larger than their al- 
ready large value due to dust). The first of these two evolved 
galaxies was best fit by a 450 Myr, 3.8 x lO'' Mq stellar popu- 
lation with a rest-frame EW of ^ 100 A. This relatively high 
EW at old age was due to dust, as the best-fit model had 
0.3 magnitudes of dust extinction at 1200 A, but this dust 
was in a geometry such that it did not attenuate Lya mea^ 
surably. This resulted in a ~ 35% increase in the Lya EW 



over a similar stellar population with no dust. The second of 
these objects had a most likely age (sec F09 for derivation of 
most likely age) of 500 Myr, a mass of 5 X 10" Mq, and 3.5 
mag of dust extinction in a very clumpy geometry, greatly 
enhancing the Lya EW. 

While these objects are very interesting, and could go 
a long way towards explaining the apparent discrepancy in 
EW distributions, the rest of the objects in the sample had 
their own surprises. They were all young and lower-mass, 
with ages from 3-15 Myr, and masses from 1.6 x 10* - 6.1 
X 10^ Mq . These ages and masses are indicative of the tradi- 
tional belief of Lya galaxies, where they are very young, and 
possibly primordial in nature. However, all of those twelve 
LAEs were best fit by stellar populations with considerably 
non-primordial amounts of dust, with A1200 = 0.5 - 4.50, 
with many galaxies showing indications of a clumpy ISM. 

By fitting their observations to stellar population syn- 
thesis models, many other recent studies have derived dust 
extinction in their Lya galaxies as well (e.g.. Chary et al. 
2005; Pirzkal et al. 2007; Lai et al. 2007; Pentericci et al. 
2008). However, these model parameters can be degenerate, 
as at rest-frame UV and optical wavelengths the reddening 
of the continuum colour could bo due to both dust and old 
stars. While observing these objects in the rest-frame optical 
with the Spitzer Space Telescope helps reduce this degener- 
acy, data at longer wavelengths could provide much tighter 
constraints. 

When UV and optical photons arc absorbed by dust, 
they are re-emitted in the far-infrared (FIR), resulting in 
a second observed peaJc in the spectral energy distribution 
(SED). Observing these galaxies in the sub- millimeter and 
millimeter (mm) regimes would detect this dust emission di- 
rectly, which could resolve this red-colour degeneracy. The 
Atacama Large Millimeter Array (ALMA) will be a uniquely 
sensitive tool to detect these high-redshift galaxies if they 
have dust, given its excellent observing conditions and large 
baseline. Thus, obtaining ALMA observations (even non- 
detections), will provide a much tighter constraint on the 
amount of dust in these objects than rest-frame UV-optical 
observations alone. We use ALMA as a measure of future ob- 
servatory capabilities of detecting dust omission from LAEs. 
If we can confirm that many Lya galaxies do in fact have 
dust, then it will help show that they are not the primitive 
galaxies that they were once thought to be. On the other 
hand, if we do not find evidence for dust, it may indicate 
that a top-heavy IMF or Population HI stars are responsi- 
ble for the strong Lya EWs. 

Where applicable, we assume Benchmark Model cos- 
mology, where fim = 0.3, JIa = 0.7 and Ho = 0.7 (c.f. Spergel 
et al. 2007). All magnitudes are listed in AB magnitudes 
(Okc & Gunn 1983). In Section 2 we discuss our sample of 
Lya galEixies which we use to make predictions. In Section 3 
we present our results, and in Section 4 we make conclusions 
about the implications of these results. 



2 METHOD 

2.1 Lya Galeixy Sample 

We examine samples of LAEs from two recently published 
studies: Finkelstein et al. 2009 and Pirzkal et al. 2007. The 
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sample from F09 consists of 14 Lya galaxies in the Great Ob- 
servatories Origins Deep Survey (GOODS) Chandra Deep 
Field - South (CDF-S) narrow-band selected to be at z w 
4.5. We will use the identifiers from F09, which are based 
on their detection filters: CHa-1 to 4 (detected in NB656, 
or Ha), CH8-1 and 2 (detected in NB665, or Ha + 80), 
and CS2-1 to 8 (detected in NB673, or [SII]). See F09 for 
further details on the data reduction and selection process. 
Seven of the 14 objects have published photometric redshifts 
putting them near z ~ 4 4.5 (MUSIC, Grazian ct al. 2006; 
FIREWORKS, Wuyts et al. 2008). As mentioned in the in- 
troduction, 12 of these objects were best-fit by young (3 - 
15 Myr), low-mass (1.6 60 x 10* M©) stellar populations 
with a good amount of dust (A1200 = 0.5 - 4.5), with the 
dust in some cases enhancing the Lya EW. Two objects 
(CHa-3 and CS2-4) were determined to contain older (450 

- 500 Myr), more evolved populations, with dust extinction 
attenuating the continuum only due to a clumpy ISM, thus 
enhancing the Lya EW. 

The sample from Pirzkal ct al. (2007; hereafter P07) 
consists of nine Lya galaxies selected in the Hubble Ultra- 
Deep Field (HUDF; Beckwith et al. 2006), identified by their 
emission lines with spectra from the GRAPES (Grism ACS 
Program for Extragalactic Science; Pirzkal et al. 2004; Mal- 
hotra et al. 2005; Xu et al. 2007; Rhoads et al. 2008) survey. 
These objects were selected based on a detectable Lya emis- 
sion line in the GRAPES grism data, and broadband ACS 
photometry consistent with a Lyman break galaxy (P07). In 
this paper, we will use the object IDs from P07: 631, 712, 
4442, 5183, 5225, 6139, 9040, 9340 and 9487. These objects 
lie in the redshift range 4.00 ^ z ^ 5.76, with Lya line fluxes 
ranging from 17-60 x 10~^* erg s"'^ cm~^ and 2' magni- 
tudes from 25.66 - 28.36. Using an exponentially decaying 
star formation history (SFH), P07 found best-fit ages of 1 

- 20 Myr and stellar masses from 7 x 10^ - 1.8 x 10® Mq. 
Three of the objects were best fit by dust, with 0.05 - 0.6 
magnitudes of extinction in the V-band (four objects were 
best-fit with dust using a single-burst SFH, and five objects 
were best-fit with dust using a two-burst SFH). 

For comparison with our LAEs, we also use the catalog 
of spectroscopically confirmed HUDF Lyman Break Galax- 
ies (LBGs) from Hathi, Malhotra & Rhoads (2008; hereafter 
HMR08; Malhotra et al. 2005; Rhoads et al. 2008). There are 
49 LBGs in the catalog, but for objects at redshifts greater 
than 4.5, we require J and H band measurements to derive 
their FIR fiuxes, thus we were restricted to using only 21 of 
the 49 LBGs. 



2.2 Dust Emission 

While we know that photons absorbed by dust are re- 
emitted in the far-infrared, the stellar population models 
used in both F08ab and P07 (Bruzual & Chariot 2003), 
while computing fluxes as red as 160 /im rest-frame, do not 
compute dust emission. In order to predict the FIR dust 
emission from the objects in our sample, we exploit a rela- 
tion between the UV spectral slope (/3; defined as ix oc A'' by 
Calzetti et al. 1994) and the infrared excess (defined as the 
ratio of FIR fiux to UV flux) detailed in Meurer et al. (1997; 
hereafter M97). This relation is shown in Figure 1 of M97, 
with the infrared excess increasing exponentially with the 
spectral slope. This relation is due to dust reddening of the 



spectral slope (/3 becoming less negative). Larger quantities 
of dust will increase the amount of UV extinction, increasing 
the amount of flux re-emitted in the FIR. 

Since we have obtained estimates of dust in our previ- 
ous studies, we could use our derived extinction values to 
estimate the FIR flux. However, not all objects are able to 
be detected in the observed IR, especially at higher redshift. 
These IR data are necessary to breaJi the aforementioned de- 
generacy of red colors due to dust or old stars. We aim to 
predict FIR fluxes independent of stellar population model- 
ing results, with the goal that our method can be applied to 
objects which only have detections in the optical or NIR. 

To use the M97 relation, we first require the UV spectral 
slope for all objects in our sample, which usually requires 
spectra. However, in M97, the authors detail a method of 
obtaining /3 from the broad-band V and R magnitudes. Un- 
fortunately, this relation is only valid for z ~ 3, but HMR08 
compute similar relations for z ~ 4 and z ~ 5 - 6. For objects 
in our sample with z ^ 4.5, we computed (3 via: 

/3 = 5.65 X (z' -z')-2 (1) 

and for objects with z > 4.5, we computed /3 via: 



' = 2.56 x{J-H)-2 



(2) 



Thirteen of the 23 objects in our sample were calculated to 

have (5 ^ -2.5, indicating that the dust in these galaxies is 
contributing significantly to their FIR fluxes. Table 1 de- 
tails the relevant parameters for each object in our sample, 
including the UV spectral slopes calculated with the above 
method. 

Once B is calculated, we use the relation from M97 to 
compute the FIR flux from the UV flux (defined in M97 
as the fiux at 2320 A). First, we used the continuum slope 
calculated in equation 1 or 2, and then from the definition 
of continuum slope (see above), we calculated the UV flux 
via: 



/a, 



/ [2320 • (1 + z)f+^ 



9200" 



(3) 



where Ia.z is the z' flux density in units of fA and 9200 A is 
the central wavelength of the z' filter. To find Ffir/Fuv, 
we fit an exponential function to the dashed line in Figure 
1 of M97. The equation of this line was: 



Ffir 1.41. a 
= 66.3 ■ e 



1.92 



(4) 



This equation accurately matched M97's Fig. 1 over the 
range of -2.5 < /J < 4. While some objects in our sam- 
ple did have (3 < -2.5, these objects will have a very small 
Ffir/Fuv ratio, and thus the emission from their small 
amount of dust will very likely not be detectable. 

In order to obtain the FIR fiux densities at a few repre- 
sentative wavelengths, we used equation 10 from M97, which 
details how the FIR fiux is composed of the flux densities at 
60 /im and 100 /xm (in Jy): 



Ffir = 1.26 X IQ-" (2.58/60 + /lOo) erg cm'^ s" 



(5) 



To independently calculate feo and fioo , we require sec- 
ond equation relating these two parameters. We used a tem- 
plate SED from Bendo et al. (2006), who computed the best 
fit modified blackbody curve to NGC 4631, an edge-on Sd 
galaxy (a reasonable local analog to LAEs). They found that 
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a simple blackbody modified by a A ^'^ term, with a dust 
temperature of 28 K, fit tlieir observed data points from ~ 
10 - 1000 ^im. 

B^iT)0C^^j^^;jj^{-) (6) 



Adopting an SED of this form and T = 28 K, we calculated 
a ratio of fioo/fso of 3.64. Using this in combination with 
Eqn. 5, we found: 



/eo 



Ffir 



1.26 X 10-11(2.58 + 3.64) 



Jy 



(7) 



3 RESULTS 

While we know the rest-frame 60 /im flux for the objects in 
our sample, we also would like to know how their FIR flux 
changes with wavelength. To do this, we scaled the Bendo 
et al. modified blackbody function up to the value of the 
60 /xm flux for each object (redshifted to the redshift of the 
object), so that we could see how their fluxes behaved at all 
wavelengths in the range of interest. We can then compare 
the flux densities of these objects at any FIR wavelength to 
the sensitivities of sub-mm instruments to see if an object 
will be detected. With the ALMA sensitivity calculatoiQ, we 
computed ALMA sensitivity curves using the default set of 
parameters for an integration time of four hours in bands 8 
(385 - 500 GHz) and 9 (602 - 720 GHz). Continuum mea- 
surements are possible with ALMA in 8 GHz instantaneous 
bandpasses, corresponding to a bandpass width of 14 yLxn 
at 730 /im (Wilson 2007). We thus chose three windows 
in which to compare our objects to the ALMA sensitivity 
curves: 450, 610 and 730 yLxn. The average sensitivity in an 
8 GHz width across these three bandpasses are plotted as 
thick black lines in Figure 1. From the relative positions of 
these sensitivities to the SED curves from our objects in this 
figure, we can see which objects will have their dust emission 
detected with ALMA. 



3.1 Lyman Alpha Galaxies 

Out of the fifteen objects in the CDF-S sample. Fig. 1 shows 
that six of them appear to contain enough dust to be de- 
tected in the observed sub-mm. These seven objects, CHa-4, 
CH8-1, CS2-2, CS2-3, CS2-4 and CS2-6, contain dust in a 
range A1200 = 1.0 - 4.5 mag. One of these, CS2-3, which 
was best-fit by A1200 = 4.5 mag, has the brightest fioo data 
point (28333 pJy) in the whole sample, further proof that 
our derivation of FIR fluxes is consistent with our model fit- 
ting. Objects CHa-1, CHa-3, CH8-2, CS2-5 and CS2-8 have 
P < -2.5, thus we did not calculate their FIR fluxes. Objects 
CHa-2, CS2-1 and CS2-7 do have measurable 60 /im fluxes 
(0.2 - 4.7 /iJy), but these are likely much too faint to be 
detectable. It is of note that CS2-4 is only detectable in the 
730 /im band. 

Three of the GRAPES LAEs which we analyzed appear 
to be bright enough in the FIR to be seen with ALMA, al- 
though only in the 730 fim band. These objects are 5183, 
6139 and 9040. In Pirzkal et al. (2007), object 5183 has a 
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Figure 1. Derived FIR SEDs versus observed wavelength of our 
objects, with 13 LAEs being shown in the plotted range. The 
different coloured lines show the LAEs, while the dotted black 
lines show the SEDs of the sample of LBGs. The curves were 
derived by scaling the computed 60 fim flux point to a modified 
blackbody curve (Bendo et al. 2006). The FIR fluxes of these 
objects were computed via the relation between the UV spectral 
slope and FIR excess from Meurer et al. (1997). The thick black 
lines show the ALMA sensitivities in a four hour integration, in 
three windows at 450, 610 and 730 fim. Out of the whole sample, 
nine LAEs (13 LBGs) would have their dust emission detected in 
a 4 hour ALMA integration. 

maximum dust extinctior0 of Av = 0.6 magnitudes (equiva- 
lent to A1200 ~ 2.4), and object 6139 is best fit by 0.4 magni- 
tudes of dust extinction. This is consistent with our results, 
which imply that the UV light absorbed by this dust is being 
re-radiated in the FIR. While object 9040 was not best-fit 
by any dust, the authors do acknowledge that its best-fit 
model was still a poor fit, although they are confident that 
it is a LAE at its measured redshift. An important point 
is that the CDF-S LAEs have a higher detection fraction 
than the GRAPES LAEs. This is hkely due to differences 
in the sample selection method. The CDF-S sample was se- 
lected via ground-based narrowband imaging (with space- 
based broadband photometry). However, the GRAPES sam- 
ple was selected via space-based grism spectroscopy over a 
much smaller area (only the HUDF vs. the whole CDF-S), 
thus probing much deeper into the faint end of the lumi- 
nosity function. As the detectability of an object scales not 
only with its color, but with its continuum brightness as 
well, fainter objects will be more difficult to detect. 



3.2 Lyman Break Galaxies 

Figure 1 also plots the derived FIR SEDs of the sample of 21 
LBGs. Dust emission from 11 of the LBGs are detectable at 
450, 610 fim and 730 fim, with two additional LBGs being 
detectable at only 730 fim. In a rough comparison between 
LBGs and LAEs, 39% of the 23 LAEs are detectable with 

^ This value of dust extinction is the maximum amount of ex- 
tinction present in this object, because the models were fit to 
this object using upper limit constraints on the rest-frame opti- 
1 Available online at: |http://www.eso.org/projects/alma/science/bin/sg ]ajifim<tgs html| 
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Table 1. 


Properties 


of LAEs 


in the Sampl(j^ 








Name 


Redshift 


/3 


fuv 


feo 


/lOO 


Dust Detectable? 








(10"^" erg cm~^ 




(^Jy) 


(Y/N) 


CHa-1 


4.40 


-3.46 


1.84 






N 


CHa-2 


4.40 


-2.50 


3.41 


0.178 


0.648 


N 


CHa-3 


4.40 


-2.90 


2.76 






N 


CHa-4 


4.40 


-1.21 


10.46 


169.1 


615.4 


Y 


CH8-1 


4.47 


-0.35 


3.97 


248.1 


903.2 


Y 


CH8-2 


4.47 


-2.67 


5.13 






N 


CS2-1 


4.53 


-2.22 


2.96 


4.733 


17.23 


N 


CS2-2 


4.53 


-0.52 


2.65 


129.7 


472.2 


Y 


CS2-3 


4.53 


1.52 


8.44 


7784 


28333 


Y 


CS2-4 


4.53 


-1.14 


0.65 


12.05 


43.87 


Y 


CS2-5 


4.53 


-4.70 


0.32 






N 


CS2-6 


4.53 


-0.10 


1.37 


125.3 


456.1 


Y 


CS2-7 


4.53 


-1.65 


0.44 


3.296 


12.00 


N 


CS2-8 


4.53 


-2.66 


0.70 






N 


631 


4.00 


-2.34 


2.28 


1.779 


6.477 


N 


712 


5.20 


N/A 


0.747 






N 


4442 


5.76 


-4.94 


0.145 






N 


5183 


4.78 


-0.75 


0.327 


11.80 


42.97 


Y 


5225 


5.42 


-3.61 


3.422 






N 


6139 


4.88 


-2.03 


4.191 


13.63 


49.61 


Y 


9040 


4.90 


-1.87 


3.093 


15.27 


55.58 


Y 


9340 


4.71 


N/A 


0.521 






N 


9487 


4.10 


-2.62 


0.822 






N 



This table lists some of the relevant properties of the LAEs in our sample. Objects with /3 
< -2.5 were not able to have their FIR excess derived via the Meurer relation, and thus we 
do not tabulate their fluxes at 60 and 100 fim (with such a steep slope, they likely would not 
have detectable dust emission). 



ALMA, compared with 62% for LBGs. While these numbers 
may be consistent given the relatively small sample sizes, it 
could be indicative of a larger difference between the types 
of objects. Traditionally, LBGs are thought to harbor more 
evolved and more massive stellar populations than those in 
LAEs. However, from this result, we now see that LBGs 
may also have a higher dust content than LAEs. This is 
confirmed by the average j3 of the two samples, which in- 
dicates that LBGs have a redder UV slope with {Plbg) ~ 
-1.62, compared with (Plae) = -2.03. A similar effect was 
seen in Vanzella et al. (2006), who found (Plbg) = -1-7 and 
(Plae) = -2.1, as well as in Pentericci et al. (2007), who 
found that LBGs were dustier than LAEs. This is interest- 
ing, given that in F09 we found that our sample of LAEs 
had much higher dust extinction rates than we would have 
thought. Thus, while LAEs appear to have some amount of 
processed ISM, the ISM of LBGs appears to be much more 
enriched, resulting in higher dust emission in the FIR. 



3.3 Robustness of Results 

We have chosen a value of fioo/feo ~ 3.64 from Bendo et 
al. (2006), which corresponds to a dust temperature of 28 
K. However, there is a large uncertainty in the FIR color. 
Larger values of fioo/feo imply a cooler ISM (thus a cooler 
dust temperature), while smaller values imply a warm ISM, 
more consistent with compact star formation (Dale et al. 
2001). Dale et al. (2001) and Dale and Helou (2002) com- 



pute model FIR SEDs assuming a large range in FIR color, 
from 0.8 < fioo/feo < 3.3. Since LAEs are expected to be 
undergoing significant bouts of star formation, and they are 
quite compact (e.g., Pirzkal et al. 2007), it is interesting to 
see how our derived FIR SEDs would change if we assume a 
warmer value of fioo/feo- The observations presented in Dale 
& Helou (2002) show that the number of objects drops off 
considerably below fioo/feo = 1, so we assumed fioo = feo, 
and redid our calculations. Figure 2 shows these results for 
CS2-2, which was one of the fainter LAES detected in all 
ALMA bands. Since warmer FIR SEDs will result in fainter 
fluxes in the ALMA observable range, this plot can give us 
an idea of the limits on our predicted fluxes. 

We flnd that although the shape of the SED has changed 
as expected, this particular object would still be detected. 
However, since the peak of the dust emission has moved to- 
ward shorter wavelengths, objects which were only detected 
in the 730 ^m band may be less likely to be detected if 
they have a warmer SED. As another illustration of our un- 
certainties, we have plotted the 1 a errors on our original 
predicted 60 and 100 fluxes, obtained from the photo- 
metric i' and z' errors. We also plot error bars derived from 
the uncertainty in the M97 relation, as described in Meurer 
et al. (1999; 0.04 dex uncertainty in log{FFiR/Fuv))- These 
errors are much smaller than the photometric errors. How- 
ever, it is worth deeper investigation into the extent which 
the photometric errors can alter our results. 

Using equations 1 - 6, we can calculate the dust emis- 
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Figure 2. Same plane as Figure 1. The blue line plots the pre- 
dicted rest-frame FIR SED for CS2-2. The literature shows that 
the temperature of the SED depends on whether the star forma- 
tion proceeds in compact or widespread environments. We thus 
show (in red) our predictions if we assumed a warmer SED, with 
T = 46K. The change in temperature does not affect the detec- 
tion of this object, although it does move the peak of the dust 
emission towards lower wavelengths. This can affect our results 
when an object is only detected in the reddest ALMA band. We 
also show the derived errors on our original SED, obtained from 
the i' and z' photometric errors. The tiny error bars inside the 
photometric errors are the errors derived from the scatter in the 
M97 relation (Meurer et al. 1999). 



Figure 3. A colour-magnitude diagram, depicting the position 
of all objects with z ^ 4.53. Filled circles represent objects which 
would be detected with ALMA at 730 /im (and possibly also at 
450 and 610 /tm), with open circles representing those objects 
which would be undetected (based on Figure 1). Blue circles rep- 
resent the CDFS LAEs, while red denotes the HUDF LAEs and 
red the HUDF LBGs. The contours show the observed fluxes at 
730 /.im for a computed grid of magnitudes and colours {il — z'). 
The black contour shows the 4 hour ALMA 730 /^m sensitivity 
limit. The contours were all computed for z = 4.33, which is the 
mean redshift of all objects in our sample at z ^ 4.53. The error 
bars on the data points are the la errors. 



sion from any object given a z' magnitude and a colour 
(i' — z for redshifts less than 4.5, and J — H for redshifts 
greater than 4.5). Figures 3 and 4 show these results at 730 
^m, with contours showing the flux levels for any colour- 
magnitude combination, and one contour speciflcally show- 
ing the ALMA 4 hour sensitivity at 730 /im. Figure 3 plots 
all objects in our whole LAE-I-LBG sample with z 5C 4.53. 
The mean redshift of this subsample is 4.33, so the flux 
curves were derived for this redshift. Figure 4 plots all ob- 
jects with z > 4.53, with the flux curves derived for the mean 
redshift of this subsample, which was 5.26. 

The different coloured circles represent the three dif- 
ferent samples (CDF-S LAEs, HUDF LAEs, and HUDF 
LBGs), with flUed circles representing objects we found that 
would be detectable with ALMA via Figure 1. All filled cir- 
cles lie above the ALMA 4 hour sensitivity contour (black 
line) showing that Figures 1 and 3 & 4 are mutually consis- 
tent, although the three detectable HUDF LAEs do lie less 
than lo" above the ALMA 4 hour sensitivity contour. Us- 
ing these figures, one could then take the colour, magnitude 
and rough redshift of any object, and determine if the dust 
emission could be detected by seeing if the points lie above 
the ALMA sensitivity contour. 

Figures 3 & 4 allow us to examine the robustness of our 
results, as we have visualized the la errors on our objects 
fluxes. We can thus see how close our objects are to the 
ALMA sensitivity limit as a function of their photometric 
errors. While 39% of the LAEs will have their dust emission 
detectable, only 17% (4/23) of them lie more than la from 
the ALMA sensitivity contour. Likewise, for LBGs, only 33% 
lie more than la from the ALMA 4 hour sensitivity contour. 




24 25 26 27 28 

z' 



Figure 4. The same as Figure 3, for all objects with z > 4.53, 
thus the colour is J - H. The contours were all computed for z = 
5.26, which is the mean redshift of all objects in our sample at z 
> 4.53. 

However, some objects which lie below the ALMA contour 
are less than la away. Thus, we can deflne a maximum per- 
centage of objects with detectable dust emission, including 
not only those above the line, but those that are below but 
within la. We flnd this to be 61% for LAEs and 81% for 
LBGs. This allows us to put an error bar on our LAE de- 
tection fraction, of 39 ± 22 % (likewise 62 Ijg % for LBGs). 

While the FIR - (3 relation was derived at low red- 
shift, both Seibert, Heckman & Meurer (2002) and Reddy 
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& Steidel (2004) find that it also appUes at high-redshift. In 
addition, Teplitz et al. (2006) found that when they com- 
pared FUV to Spitzer observations, the resuhs were in gen- 
eral agreement with the FIR - /3 relation. Chapman et al. 
(2005) have investigated this relation using a sample of sub- 
millimeter selected galaxies at z ~ 2. The authors conclude 
that the bolometric luminosity of sub-millimeter galaxies 
are underestimated by 2 orders of magnitude when derived 
from their UV properties. This indicates that our derived 
FIR fluxes may be lower limits, meaning that if anything we 
could be underestimating the detection fraction. 

However, Reddy et al. (2006) found that the M97 re- 
lation may not apply to all star- forming galaxies at z ~ 2. 
They find that galaxies with derived ages < 100 Myr ex- 
hibit a smaller FIR excess than that predicted by the M97 
relation, indicating that these young galaxies have redder 
UV colors for a given amount of dust extinction that older 
galaxies. This is of concern as most of our LAEs have very 
young derived stellar population ages. In addition. Baker 
et al. (2006) find that the M97 relation also over-predicts 
the FIR flux for the gravitationally lensed LBG cB58. How- 
ever, there are differences between these objects and LAEs, 
as they are at a lower redshift, and they are continuum se- 
lected. We thus cannot know if the same age dependence will 
apply to LAEs until we obtain a statistical sample of rest- 
frame FIR detections of LAEs. In addition, the Reddy et al. 
flgure which illustrates the age discrepancy does not have 
error estimates on their young object points, thus it is hard 
to discern how signiflcant this effect is. Baker et al. do show 
error bars, and these errors imply that their observations 
are only off the M97 relation by ~ 1 cr. 

If young objects do have systematically lower 
Ffih/Fuv ratios, it would be worrying, as it could imply 
that we are over-predicting the FIR fluxes of our LAEs by 
as much as an order of magnitude. It is thus prudent to 
apply our results to any relevant observations to see if our 
FIR flux predictions are accurate. Although very few galax- 
ies with Lya emission have been detected in their rest-frame 
FIR, Capak et al. (2008) have recently published a study of 
a spectroscopically conflrmed millimeter starburst galaxy, 
which does exhibit Lya emission. They publish the observed 
SED of this galaxy from X-rays to radio wavelengths. Using 
the published optical magnitudes, we can compute /3, and 
predict the FIR flux. Doing this, we calculate a predicted 
554 (rest 100) /im flux of 2637tlatl ^Jy- However, the clos- 
est point in the observed SED is at 1.1mm, at 4800 ± 1500 
/iJy. The flux of our derived SED at 1.1mm is 20541^^^^ 
^Jy, thus they are consistent within 1 a. These observations 
give evidence that our predictions do match observations. 
We show the SED from Capak et al. along with our predic- 
tions for this object in Figure 5. Additionally, we also show 
the V, i', z' and 3.6 fim fluxes for CH8-1, which was one of 
our redder LAEs. These points have been scaled to match 
the i' flux of the Capak et al. object. We can see that the 
i' - z' color is somewhat redder than the Capak et al. ob- 
ject (0.3 versus 0.15), which indicates that CH8-1 is dustier, 
explaining its higher expected FIR flux. However, as it is 
intrinsically fainter (the scaling factor was ~ 17), it would 
be much harder to detect with current instrumentation. 

Although we have assumed a reasonable SED shape, 
it is possible that the FIR SED could be much different. 
Houck et al. (2004) discuss the low metallicity blue compact 
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Figure 5. In this figure we plot the SED from Capak et al. (2008) 
for a millimeter detected Lya emitting starburst galaxy, spectro- 
scopically confirmed to lie at z=4.547. The predicted FIR fluxes 
match up well with the observed SED, providing further evidence 
that the M97 relation is applicable at high redshifts. The blue 
points denote the V, i', z' and 3.6 fira fluxes from CH8-1. These 
points have been scaled to match the i' flux of the Capak et al. 
object. We can see that the i' - z' color is redder than the Ca- 
pak et al. object, which indicates that CH8-1 is dustier, which we 
can see as its (scaled) predicted FIR flux is brighter than that of 
the Capak et al. object. However, this object is fainter by about 
a factor of ~ 17, thus it would be much harder to detect with 
current instrumentation. 



dwarf galaxy SBS 0335-052, which has a FIR spectrum much 
different than typical starburst galaxies, peaking at ~ 28 
Hm, indicating much warmer dust, possibly due to the lower 
metallicity. High-redshift LAEs may be expected to have 
a very low metallicity, thus it is worth mentioning that if 
any of our objects had this type of SED, we could be vastly 
overestimating their FIR flux. We have thus investigated the 
dependence of our results on these derived FIR fluxes. To 
do this, we reduced our FIR flux to be 50% of its derived 
value, and computed the ALMA detection fraction of both 
types of objects. We calculated a detection fraction of 26% 
for LAEs and 62% for LBGs (57% and 76% for LAEs and 
LBGs respectively within la of the ALMA contour). Thus, 
while we believe that the M97 relation is applicable to our 
objects, we find that our conclusions are robust to a factor 
of two variation in the FIR flux. 



4 CONCLUSIONS 

Directly observing dust emission from Lya galaxies will pro- 
vide undeniable proof that these earliest of galaxies have 
undergone some enrichment of their ISM. While these ob- 
servations are currently difficult, we have been able to pre- 
dict whether a sample of 23 LAEs would be detected with 
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a future sub-mm observatory such as ALMA based on their 
optical and NIR fluxes. Exploiting a relation from Meurer 
et al. (1997), which relates the far-infrared excess to the 
UV continuum slope, we have derived the FIR fluxes from 
these LAEs, which will be primarily due to dust emission. 
By comparing these fluxes to the sensitivity of ALMA at 
450, 610 and 730 /im, we have determined that 39% (9/23) 
of our sample of LAEs should be detectable based on their 
dust emission. 

Wc perform the same analysis on a sample of 21 Lyman 
Break Galax;ies, finding that 62% (13/21) of the LBGs are 
detectable with ALMA, ~ 60% higher than the detection 
fraction of LAEs. This may be a function of our small sam- 
ple size, or it may be the result of a true dicfiotomy in the 
dust contents of LAEs and LBGs, providing further proof 
that LBGs are a more evolved class of high-z galaxy. In or- 
der to generalize our results, we computed the FIR fluxes for 
a wide range of continuum magnitudes and colours, compar- 
ing these fluxes to the ALMA 4 hour sensitivity limit. Thus, 
knowing the continuum magnitude and colour of any ob- 
ject, one can then predict whether its dust emission will be 
detectable. 

While stellar population modeling gives some insight on 
the amount of dust in high-redshift galaxies, future sub-mm 
instruments will allow us to directly observe the dust emis- 
sion from these objects, providing a more complete census 
of dust. From our analysis of LAEs and LBGs in the HUDF, 
we calculate a number density of ALMA detectable LAEs 
and LBGs of 0.26 arcmin"^, and 1.13 arcmin"'^ respectively. 
Thus, even in a small survey area, one could expect to detect 
a number of high-redshift objects with ALMA in a reason- 
able amount of time. In our previous analysis of LAEs, the 
addition of the dumpiness parameter allows strong Lya flux 
to an arbitrarily old and dusty galgixy, thus these ALMA 
observations will help to constrain the amount of dust, and 
thus the level of dumpiness. 

These results will have considerable impact on several 
areas of galaxy evolution, specifically on the star formation 
rate (e.g., Giavalisco et al. 1996), which, being derived from 
UV observations, is critically dependent on accurate mea- 
surements of the dust content. Combining ALMA with UV 
observations, we will have a complete map of star formation 
in high-redshift galaxies (i.e. Gao 2008). Studying the sub- 
millimeter emission of these objects will also allow calcula- 
tion of the total dust mass, which is something one cannot 
calculate from UV-optical data alone. Along with learning 
about the total amount of dust in these objects, ALMA will 
give us the opportunity to investigate the distribution of 
dust due to its excellent resolution. At 730 GHz, ALMA can 
achieve resolutions down to 0."01 pix~^. Combining ALMA 
observations with high resolution HST data will allow us 
to investigate dust on local scales at high-redshift for the 
first time. Strong sub-millimeter detections combined with 
strong Lya emission could indicate that dust enhancement 
of the Lya EW is widespread. 
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